SEARCH FOR SOLITONS IN TWO-HIGGS 
EXTENSIONS OF THE STANDARD MODEL * 



T. N. TOMARAS f 

Department of Physics, University of Crete, and Research Center of Crete 
714 09 Heraklion, Crete, Greece 
E-mail: tomaras@plato.iesl.forth.gr 



in 
in 

(N 



in 



ABSTRACT 

We report on the status of our search for quasi-topological solitons of various 
dimensions in realistic field theoretical models of condensed matter and of ele- 
mentary particle physics. 



1. Introduction 



The search for stable lumps in the Weinberg-Salam model has a long hi 
It has revealed a richpgtaucture of classical ^iu4ipa£-|including the sphaleron^ 
deformed sphaleronsH'BQ'B and vortex stringsa'^tiii3li2l. Such solutions could play a 
role in understanding (B+L)-violation and structure formation in the early universe, 
O '. but they are all classically-unstable or/and extended. They have therefore no direct 

Q\ ■ present-day manifestation, contrary to long-lived particles whose relic density could 

at least in principle be detected. 

The existence of particle-like excitations has, o^the^ther hand, been argued for in 



the context of a strongly- interacting Higgs sectoBtlHllZflla The advocated particles 
can be thought of as technibaryons of an underlying technicolor model. They are 
described in bosonic language by winding solitons of an effective non-renormalizahle 
lagrangian for the pseudo-goldstone-boson (or technipion) field, much like skyrmionsi^ 
of the effective chiral lagrangian of QCD. This is of course a phenomenological 
description, since the properties of such hypothetical particles cannot be calculated 
reliably within a semi-classical expansion. Furthermore, in view of the difficulties 
^ ■ facing technicolor models, the possibility of a strongly-interacting Higgs sector is not 

H \ theoretically appealing. 

It would be clearly more interesting if classically-stable winding excitations could 
arise in a weakly-coupled scalar sector. To be more precise let us decompose the 
Higgs-doublet field into a real (positive) magnitude and a group-phase: $ = FU, 
and consider static configurations with U(x) wrapping iV times around the SU(2) 
manifold. These are potentially unstable for at least three distinct reasons: (a) 
because N is not conserved whenever the magnitude F goes through zero; (b) because 
N is not gauge-invariant and can, in particular, be non-vanishing even in a vacuum 
state; and CcJ because scalar-field configurations can loose their energy by shrinking 
to zero sizeEl. We refer to these for short as the radial, gauge and scale instabilities. 
They can be eliminated formally by (a) taking the physical-Higgs mass van — * 00, (b) 
decoupling the electroweak gauge fields, and (c) adding appropriate higher-derivative 
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terms to the action. The question is whether classical stability can be maintained 
while relaxing the above conditions. This has been investigated numerically in the 
minimala&se of one doublet: although one may indeed relax both the weak gauge 
coupling^ and the Higgs mass" up to some finite critical values, stability cannot 
apparently be achieved without the non-renormalizable higher-derivative terms in the 
action. On the other hand, as we have demonstrated recently, metastable winding 
solitons do arise in renormalizable models in twoeHl and threeE3 space-time dimensions. 
The way this happens is the content of section 2. It is we believe instructive and could 
guide the search for lower dimensional structures in many physical systems as well 
as for similar semi-classical solitons in four dimensions. Section 3 contains a brief 
presentation of the status of our search for particle-like solutions in the context of 
two-Higgs extensions of the standard model. We close with some comments in the 
discussion section. 

2. Quasi-topological solitons in lower dimensions 

The simplest context in which the radial instability is an issue is a two-dimensional 
model of a complex scalar field with mexican- hat potential: V = ±A($*$ - v 2 ) 2 . The 
simplest way to describe in detail the winding solitons is to take space to be periodic 
with period L a . The condition for classical stability can in this case be derived 
analytically and readsS : m^L > , where m# = y/2Xv. The classically-relevant 
parameter is thus the radial-Higgs mass in units of the soliton size. This follows 
also by comparing the loss in potential energy to the gain in gradient energy when 
trying to undo the winding by reducing the magnitude of the scalar. Note that the 
loop-expansion parameter XL , can be taken to zero independently so as to reach a 
semiclassical limit. 

The above winding solitons become unstable classically if we gauge the U(l) sym- 
metry of the model. The gauge instability is in fact more severe than in four dimen- 
sions, because no energetic barrier opposes the turning-on of a static space-like gauge 
field, which is necessary to reach a winding-vacuum state. The minimal Abelian- 
Higgs model has thus only unstable (sphaleron) solutions^ b . The situation changes, 
however, drastically if there are more than one complex scalars. The gauge-invariant 
relative phases of any two of them cannot in this case wind around non-triiually m 
a vacuum state. An explicit analysis of this extended abelian-Higgs modeled shows 
that winding solitons persist down to scalar masses close to the inverse soliton size: 
gauging and the extra Higgs enhance the stability region found in the global model. 

The scale instability becomes an issue for the first time in three space-time dimen- 
sions. To be more precise we consider a real-triplet scalar field Q a {x) {a> = 1, 2, 3) with 
mexican-hat potential : V = |A($ a $ a — v 2 ) 2 . The limit A — > oo corresponds to the 
0(3) non-linear a-model. This is known to possess winding solitons, characterized 
by non-trivial mappings of the two-sphere onto itself, and having arbitrary sizec^l. 
For finite A on the other handy, or in the presence of a symmetry-breaking poten- 
tial, Derrick's scaling argumentE3 shows that these solitons are unstable to shruaking. 
One can of course again invoke higher-derivative terms to stabilize the scale£3. The 

"Alternatively we may add a mass term: 8V = — /i 2 ui?e($), that lifts the vacuum degeneracy. Stahlfi 
winding excitations, which reduce to the sine-Gordon solitons in the A — > oo limit, can be showncil 
numerically to exist for Xv 2 / n > 18.8 . 

h It was claimed erroneously in [22] that it has no static solutions whatsoever. This is only correct 
in the A — > oo limit. 



same result is however inthis case achieved by a massive U(l) gauge field with only 
renormalizable couplings^. This can be established by perturbing around the 0(3) 
non-linear a- model limit, or else by solving numerically the equations of motion. 

What do these lower- dimensional solitons teach us? First, they suggest by anal- 
ogy that classically-stable winding solitons may exist in a weakly-coupled two-Higgs 
extension of the standard model. The status of our search for such localized solutions 
is the content of the following section. 

Second, they are interesting in their own right, since they correspond to a new 
class of wall and string defects in realistic four-dimensional models of particle and 
condensed- mattpi r«hysics . One such example are the metastable membranes, dis- 
covered recentlycS'E] in two or more Higgs-doublet extensions of the standard model. 
These are static wall-type solutions, non-topological but classically stable in a wide 
region of the Higgs sector parameter space compatible with perturbative unitarity and 
with present phenomenological bounds. They are embeddings of the above discussed^ 
solutions of the 2d Abelian-Higgs model with two or more complex scalars. They are 
characterized by the non-trivial winding of the relative U(l) phase of the neutral 
components of any two Higgs doublets of the extended standard model in the direc- 
tion x normal to the wall. They have no electromagnetic coupling, their size is a few 
times the inverse of the mass mj of the CP-odd scalar A (in the standard notation 
of the generic two-Higgs modeled) , while their energy per unit area is in terms of the 
A and W masses and the fine structure constant of order m^m^/a. Assuming a 
m A ~ 50GeV i.e. not much larger than the present experimental lower boun d§, the 
mass of the wall is of O(10 10 gr /cm 2 ). Thus, a single wall crossing the entire Universe 
today would by far overdose it and can be excluded. Smaller membranes though, 
which either collapsed or were torn apart by quantum tunneling may have acted as 
seeds for the formation of galaxies. In fact, the mass of a typical galaxy is comparable 
to that of a membrane a few light years in size. In any case, the production and decay 
rates of these objects and their cosmological role has to be studied in detail, since 
most likely they are going to lead-to firm constaints on the Higgs mass spectrum of 
extensions of the standard modeH 

As another example, we would like to mention the static wall-type metastable 
solutions of the easy-plane ferromagnetic continuum in the presence of an external 
in-plane magnetic field h c . The dynamics of the system is described by the Landau- 
Lifchitz equations with energy density given in terms of the unit magnetization vector 
n a , a= 1,2, 3, by: £ = \(din a ) 2 + gn% 2 + h{l — ni). The anisotropy constant g and the 
external field h are both taken to be positive. The system has the unique semiclassical 
ground state: {rii = l,n 2 = = n 3 ) and no topologically stable domain walls. 
Nevertheless, it posesses a variety of wall-type static finite energy (per unit area) 
solutions. A particularly interesting one is {n\ = cosQ(x),n2 = sinO(x),n 3 = 0), 
with 6(x) satisfying Q"(x) = hsinO(x) and the boundary conditions G(— oo) = 

c The search for finite-energy soliton solutions in ferromagnetic and antiferromagnetic systems and 
the study of their stability, reported briefly in this paragraph, was carried out in collaboration with 
Prof. N. Papanicolaou. I would also like to thank Dr. P. Tinyakov for helpful discussions about 
sphaleron solutions in relativistic anisotropic non-linear a— models in 1+1-dimensions. 



0, 0(+oo) = 2tt. It can be verified that this solution is dynamically stable for g/h > 3. 

3. Localized solutions in 3+1 dimensions d 

We now turn our attention to the possibility of classically stable localized particle- 
like winding solitons in 3+1 dimensions. In the context of a two-Higgs extension of 
the standard model these hypothetical solitons would: (b) be characterized by the 
non-trivial winding of the relative phase of the two doublets, and thus be immune to 
the gauge mode of decay; (c) have a scale stabilized by electroweak magnetic fields and 
hence of order 1/mw', and (a) hopefully stay stable for Higgs masses near mw and thus 
compatible with perturbative unitaKvty e . Mathematically the situation is the same as 
in the hidden-gauge-boson models"'^ of strong and electroweak interactions, except 
that the role of the hidden gauge bosons, is, -here played by W and Z themselves. 

We carried out a numerical searcht3a'EH for such stable spherically symmetric 
particle-like soliton solutions in a simplified version of the two-Higgs SU(2j x U(l) 
model. The presumably inessential simplification consists of taking (a) the U(l) gauge 
coupling g' = as required by spherical symmetry and (b) the Higgs potential to be 
the sum of two "mexican hats" one for each Higgs doublet. The resulting action is: 

S = i / *x(-±Tr{W ia ,W' a ')+ E tP l >H 1 jt{iy>H I ) + E ^(HjHj-g^r) (1) 

9 J Z 1=1,2 1=1,2 9 

with D^Hj = (dp + W^Hj, 1 — 1,2 the covariant derivative on the Higgs doublets. 
Wfi = ^ , r a H / ° and W^ v is the SU(2) field strength. We used the spherically symmetric 
ansatz 

W = a Tini/2i , 
Wi = [(a - VjtijkTjTik/r + (3(8ij - riiU^Tj/r + a 1 n i n j T j }/2i 
Hi = (iii + wmTi)^ 

where rii is the unit vector in the direction of x, r» are the three Pauli matrices, ^ 
is a constant unit doublet and do, ai, a, f3, fij, i/j are functions of r to be determined 
by the energy extremization. It is convenient to choose the gauge ao = 0. One then 
solves Gauss' constraint to obtain a\ = 00. 

We rescale distances and fields by the appropriate powers of mw = gJ (v\ + f|)/2 

and use the notation 0/ = /// + ivi = Fie 1 ® 1 to write for the energy functional: 

E = j dr{^f(a" + P") + r 2 (Ff + if + F«Q» + F^) 

+ ^(« 2 + P 2 ~ I) 2 + \(* 2 + P 2 + 1)(^ 2 + F*) 

d The work in this section was done in collaboration with Dr. P. Tinyakov and is described in detail 
in reference [31]. 

e Though admittedly premature, some other physical properties of such would-be particles are fun 
to contemplate: being classically stable they could easily have cosmological life times. They would 
have a mass in the ~ 10 TeV region, zero charge and dipole moments in their ground state, and 
geometrical interaction cross sections of order Assuming maximum production at the elec- 

troweak phase transition, a rough estimate of their present abundance shows that they could be 
candidates for cold dark matter in the universe. 



-a[F 2 008(20!) + Fl cos(20 2 )] - /3[F? sin(26i) + F 2 sin(20 2 )] 



+ fr 2 (F 1 2 -l) 2 + ^r 2 (F|-t 2 ) 2 }, (2) 

where t = v 2 /vi and k i = m Hl /m w . The radial coordinate in particular in the above 

expression is measured in units of m^. 

We require finiteness of the energy of the solution and use the global symmetries 
of the energy functional^ to specify the boundary conditions at oo and write them in 
the form: 

a->l, /3->0, Fi->1, F 2 -*t, a ^O 
Correspondingly, at r = energy finiteness implies 

a 2 + /? 2 -l. 

while the field equations obtained by varrying (3) with respect to 0/, together with the 
requirement of smoothness of the solution lead to the dynamical conditions: 20 a — ip = 

mod 27r , where ip is the phase of the field x = a + 0- Thus all solutions satisfy 

01 — ©2 — at the origin. 

We have looked for solutions in a wide region of the {t, k±, k 2 } parameter space. 
We refer the reader to [31] for the detailed presentation of the numerical method we 
used and of the entire solution-zoo we have obtained so far. In the present note, I 
will describe briefly some of the solutions we found for various values of k on the 
line {t — 1,k± — K2 = k}, together with their main characteristics. As we varry k 
the energy landscape changes smoothly and consequently, one generically expects a 
continuous change in the number and properties of its extrema and saddles. Figure 1 
below captures these changes in part of the energy landscape and shows some of the 
solutions of the model at hand. 

The graph should be thought of as an "artist's conception" of a part of the tree of 
solutions with N = and \N\ — 1 which has emmerged so far in the two-Higgs model 
under study. Each branch of the "tree" corresponds to a particular type of solution. 
As we move along a branch by varrying k, continuous quantities like the detailed form 
of the solution, its energy as well as the value (s) of the imaginary eigen-frequencies 
of the small oscillations around it change smoothly, while integer-valued ones such 
as the index N, the number of unstable negative curvature modes (shown on the 
branch) and the number of nodes of the fields, remain invariant. We have suppressed 
the branches corresponding to the CP-conjugates of the solutions shown. They may 
be imagined as the mirror images with respect to the k— axis of the branches above 
it. 

For k between 0.0 and 5.5 the only non-trivial solution is the sphaleron embedded 
in the two-Higgs model; it has ©i = 2 = 0, (3 = and one negative mode. For k 
about 5.5 we reach the bifurcation point B , the sphaleron develops a second negative 
mode and at the same time two new solutions, namely the branch A 1 and its CP 
conjugate appear with N = 1 and N — — 1 respectively. The phases ©i and 2 and 
the complex field x °f th e solution corresponding to the branch A 1 have the behaviour 
we call of "type- A" and show in figure 2. 

^The energy functional is invariant under the simultaneous rotation of the complex fields 
0i,02 and a + i[3 by an angle u,u) and 2uj respectively, as well as under independent rotations 
of 4>i by 7r. 



We have checked that the branch Ai extends without a new bifurcation up to 
K = 200. As for the sphaleron, once we get to k ~ 12.0 we encounter the bifurcation 
Bi. The sphaleron aquires a third unstable mode and two new solutions (branch S2 
and its CP-conjugate) emerge. They have N = and phase behaviour of "type-S" 
shown above. And so on for the rest of the tree. 

Incidentally, using the number of negative modes A neg mo des of the solutions shown 

in figure 1 one may explicitly verify the sum rule: Dsoiutions(~ ) 7Vncg - modcs = constant 
independent of k, for the whole range of values we considered. 

4. Discussion 

It is clear that the two-Higgs system discussed here posesses a variety of unstable 
solutions, whose richness increases with the Higgs masses. These, like the ordinary 
sphaleron, are in principle interesting in their own right especially in connection 
with the baryon number generation at the electroweak phase transition. There are 
important differences^ in the values of the energies and of the negative modes of 
the above solutions compared to analogous ones of the minimal standard model. For 
instance, the value of the negative curvature along the A\ branch is about half of 
what it is in the corresponding least unstable deformed sphaleron of the one-Higgs 
model. These differences affect directly the prediction for the baryon number in the 
Universe, and at the same time they offer support to our general arguments that an 
extended Higgs sector improves the stability of the winding solitons under discussion. 

Despite of this indirect evidence though, we have not as yet been able to find a 
stable soliton. Such a solution might for example arise at a bifurcation point like X 
of figure 1 and either merge with another branch of the tree at some new bifurcation 
or continue, as is the case with A\, up to k — 00. But since our numerical method 
requires an initial guess for the solution which has to lie inside the basin of attraction 
of the corresponding local minimum of the energy functional, our failure to find a 
stable solution might just be due to the bad starting configurations we tried so far. 
Further theoretical analysis is required to guide our numerical search. 
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